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Abstract: The coUisional drift wave instability is reexamined taking into account the 
ion response in the direction parallel to the magnetic field lines, which appears due to fric- 
tion with electrons and which can not be omitted in view of the momentum conservation. 
A modified instability threshold is obtained. In plasmas with dominant electron collisions 
with neutrals, the instability threshold is shifted towards higher frequencies, compared to 
the case of dominant electron collisions with ions. The difference between the two cases 
vanishes when the ion sound response is negligible, i.e., when the instability threshold 
disappears, and both ions and neutrals react to the electron friction in the same manner. 



PACS Numbers: 52.35.Kt, 52.30.Ex, 52.20.Fs 



I. Within the standard fluid theory, the drift wave can be excited by electron coUisions. 
In this case, the usual relation between the electron perturbed density and the perturbed 
potential, ni/no — e^i/^Te, becomes modifled due to the presence of the coUisional term, 
so that the potential lags behind the density.^'^ The effect appears regardless whether the 
electrons collide, in their predominant motion along the magnetic fleld vector, with ions or 
with neutrals.'^'"' Compared with the kinetic instability (which is due to the inverse electron 
Landau damping effect, that appears because the mode frequency is slightly below the 
diamagnetic frequency), the coUisional instability is dominant^ provided that the electron 
parallel mean-free path is smaller than the parallel wavelength. Hence the interest in the 
drift modes with very large parallel wave-lengths and relatively short perpendicular wave- 
lengths. In view of the early theoretical prediction^'^ and experimental veriflcation,^ the 
amount of literature dealing with the drift wave is enormous. More recently, coUisional 
drift waves have been studied experimentally in cylindric plasma in Refs. 9-11. Details 
on the experimental investigation of global coherent structures associated with the drift 
mode in simple magnetized torus can be found in recent Refs. 12-14. The presence of 
charged grains on the drift wave in a cylindric configuration has been studied in Ref. 15, 
while the role of the drift wave in the transport phenomena may be found in the most 
recent Refs. 16, 17. 

Some convenient approximations that are made in the derivations include the limit 
in which the complete dynamics of the heavier particles (i.e. ions and neutrals) in the 
direction along the magnetic fleld is negligible, and the limit when electrons can be treated 
as inertia-less 



In the presence of collisions, the friction force term in the electron parallel momentum 
equation is usually written in the form —■me'rieiyej'Ve, where j — i, n. As a result one obtains 
a standard phase shift in the electron Boltzmann distribution ni/riQ — {e(f)i/ KTe){l — iS), 
that is responsible for the mode growth. 

However, the conservation of momentum implies that the friction term in the electron 
momentum equation should read —merieiyejive — Vj) even when conventional criteria for 
a negligible parallel dynamics of heavier particles are fulflUed, so that the corresponding 
momentum component of the heavier species includes the friction term —mjnjVje{vj — Ve)- 
Below, we perform derivations with such a full friction force term for some simple cases in 
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order to demonstrate the differences introduced by tfiis friction-induced response of the 
heavier species. 



II. We first discuss a plasma with dominant collisions between the charged particles. 
Note that this case may also include a rather weakly ionized plasma with uq <S n„ (the 
index here and below denotes the electron or ion equilibrium quantities) . This is because 
of the much larger cross section for collisions between charged particles. To have dominant 
collisions with protons in a plasma containing electrons, protons and neutral atoms, the 
electron number density should satisfy the condition 

rin (87r)V2e4Le, ' ^ ' 

Here we used the standard notation, and Lgj is the Coulomb logarithm. Taking as an 
example Tg = 10^ K, which gives^^'^'^ aen = 2.5-10-^'^ m^ (here Lei = 6 for no = 10^^ m^^), 
it turns out that the electron-ion collisions are more frequent than the electron-neutral 
collisions provided that no/n„ > 0.009. For the given temperature this is close to well 
known estimate^'' showing that, in terms of electron collisions, an ion is equivalent to 
3.4 ■ 10^(300/Te)2 ~ 300 neutral atoms. 

We use the continuity equation for electrons and ions placed in an external magnetic 
field Bq — BqCz 

+ W±{njV±j) + Vz{njVzj) =0, j = e, i, (3) 
where the linearized perpendicular velocities of electrons and ions are given by 

1 - ^ '"re^ V_Lnel 2 ^ / , a\ 

v±ei = s-e^ X V±0i - —e, X , = KTe/nie, (4) 

-Do i le 

Here, the effects of collisions on the perpendicular electron dynamics is neglected in con- 
trast to the parallel one, which is justified^^ as long as klVil/ {kyuld > 1- 
The electron parallel velocity is determined from 

(901 dni 

= 6^0-^ - '^■^e-^ - menoUeAVezl " Vizi), (6) 



and the ion velocity from 

dvizi 



a '^ie{Vizl-Vezl)- (7) 

at rui az 
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(9) 



III. In what follows, we first present the 'standard' derivation in which the ion parallel 

velocity is only due to the parallel component of the perturbed electric field, while the 

friction induced term in (7) is omitted with the usual excuse of the huge difference in 

mass between the two species. In this case, the perturbed ion number density is described 
byi'2 

^-f^ + ^-^>^)^. (8) 

no \U! UJ^ "J Klf. 

Using Eqs. (3), (6) the perturbed electron number density can be written as 

no 00 + ik^Dz nTe ' 

'^Te , , kyKTg nQ 
Uei eBo no 

Here, we have taken Vno = e^dno/dx, the perturbations are assumed to be of the form 
~ f {x) exp{—iu!t + ikyU + ik^z), and we work in the frame of a local approximation. 
The full coUisional term, with perpendicular electron collisions included, should read 
i{k1v'^^/vei + plkyUei): though the second term here is negligible in the limit discussed 
earlier. Here, Pe — Vtc/^c is the electron gyro-radius. 

Note that klc^/uj in Eq. (9), that appears due to the term i^eiVizi in Eq. (6), is usually 
omitted in standard derivations [2]. The quasi- neutrality yields the dispersion equation 

^*,ky, ^^^^^ _ u, + kyju + ikiDz , . 

- + ^ - I'yPs - ^ + ,klDz ■ ^ 

We now take the limit 

< klDz, (11) 

and assume that uj and o;* are of the same order. Used for convenience^, the condition (11) 
is in fact not always easily satisfied. Physically it describes^ the condition of isothermal 
electrons along the field lines that has been assumed. It can be rewritten as {ijj/kz)/vTe ^ 
kz^Tel'^ei- The right-hand side gives the ratio of the electron mean free path VT^I'^ei and 
the parallel wavelength, that in fact must be much less than unity in order to remain within 
a proper fluid theory (i.e., for collisions being able to maintain Maxwellian distribution). 
A few comments here are noteworthy. 

i) If the parallel ion dynamics is completely neglected^ we can write the dispersion equation 
as follows: A(a;,/Cj^) = Ar(a;, /Cj,)-|-iAj(a;, /c^^) ~ lS.T.{LOr.ik^+iLOidiS.r{uj ^k^ j duj\i^=^^ + 
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iAi{u!r,ky) — 0, where |Aj| <^ jA^] and where the subscripts r,i denote 'real' and 
'imaginary'. Setting the real and imaginary parts of this dispersion equation equal 
to zero yields 

u;rc^u;,/{l + klpl), uji ~ -Ai{u;r,ky)/{dAr/du;)^=^^ ^ ^^^^ + fcV)^ " ^^^^ 

The increment is proportional to Ugi, which implies that, with the fixed ion back- 
ground (in the parallel direction), the more coUisions the larger the growth of 
the wave is. Regarding the dependence on kyPs, the increment is maximal at 
kyPs = 

ii) With the ion sound response taken into account (assuming that k^cl is same order or 
not much larger than cu^ and a;*a;), the dispersion equation becomes 



u'il + klpl) - - klcl + ^ - u.) - klcl] = 0. (13) 

Here, in the expansion on the right-hand side in Eq. (10), the sound contribution 
yields a real and an imaginary term. The former, yielding a negligible correction 
to the real frequency, is neglected, and we have only kept the imaginary term. The 
real part of the frequency is determined from 

ArioJr, ky) =ujI{1 + klpl) - u^u - klc] = 0. (14) 

This is used in the imaginary part, yielding 



klD,uj'',{l + klpl) + klc 



iii) We note that in fact in the standard literature [2] the sound contribution k^cl/uj on 
the right-hand side in Eq. (10) is neglected. It is easily seen that the origin of this 
term is the term merio^'eii'jzi in Eq. (6). In this case instead of Eq. (15) we obtain 



, ,2 /,,2i-2^2 _ jL2 2 

klD, 00^(1 + ky,) + k^,c. 



,2 



(16) 



Hence, there appears to be a threshold for the instability. 
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IV. Of course, a self-consistent analysis should include the full friction force effect in 
both, i.e., the ion and electron, parallel equations. This is simply due to the fact that the 
two forces are necessarily equal by magnitude. Hence, we keep Eqs. (6, 7) as they are, 
with the complete friction terms. 

Combining the two parallel momentum equations, and using the conservation of mo- 
mentum rriiUie = niei^ei, yields Vi^i = kzC^gUi/ {ujuq). Instead of Eq. (8), we now have 



A procedure similar as earlier, now yields the electron number density 

(18) 



ni a;* + ik'^D^ e0i 



no cu — k1c1/uj -\- ik^Dz kT^. ' 

instead of Eq. (9). Compare these two equations with Eqs. (8) and (9). Notice that 
in both Eqs. (17) and (18), the ion parallel response appears in a completely different 
manner. 

Within the same approximations as earlier (i.e., a;^,a;*a;, /c^c^ <C k^D^), the dispersion 
equation that we now have is: 

u;\l + klpl) - a;.^ - + ' ' ^^^'^ = 0. (19) 

The real frequency is the same as in the earlier Eq. (14). Hence, using Eq. (14) the 
imaginary part of the frequency now becomes 

'^^^ klDz ^^{l + klpD + klcl klDz ul{l + klpl) + klcl- ^ ^ 

We remark the obvious difference in the instability threshold in the expression (16), and 
the correct expression (20), and we stress the absence of the threshold in Eq. (15). 

V. The ion response to the friction is usually neglected on the basis of the huge difference 
in mass of the different species. However, this difference in mass may be compensated by 

the frequent electron collisions with ions, so that sooner or later the ions start to move in 
the parallel direction due to the electron drag. To get a feeling on the effects of collisions 
and the corresponding time scales, we may discuss the following two separate cases. 

a) Assuming that the condition (1) is satisfied, the ions respond in the parallel direction 
only through the friction. The electron velocity Vq is assumed nearly constant 
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by the parallel electric field of the wave. Prom (7), assuming that the ions are 
initially at rest, the ion velocity, normalized to Vq, becomes 1 — exp(— i/jet). Taking 
no = 10^^ m-^ Te = 10^ K, = 2 • 10^ K, we have ^ 7 ■ 10^ Hz. A simple 
plot of the ion velocity reveals that it becomes close to 1 already after about 0.0007 
seconds. 

b) Taking another extreme case where electrons initially, due to any external reason 
acquire a velocity Veo — Vq, without any additional force, and where Vio — 0. In 
this case the electron velocity is not kept constant, the interaction of the two fiuids 
yields the evolution of the two velocities: 



Here, the electron and ion velocities monotonously change in time towards the com- 
mon velocity (the first term on the right-hand side) Vc — Vomg/mj <^ Vq which, 
for the same parameters as above, is achieved within the time interval shorter than 

10~^ sec. The characteristic time for the velocity relaxation is ~ l/(t'ei + i^ie)- This 
is presented in Fig. 1. 

A real physical situation, as in the case of the drift wave discussed earlier, is expected 
to be somewhere in between the two extremes presented above. Hence, in spite of a huge 
mass difference, the collisions (friction) will force ions to move along the magnetic field 
lines, and, due to the same reason, the electron velocity amplitude associated with the 
drift wave is expected to be considerably smaller. 

VI. We shall check now the case of a plasma with dominant electron collisions with 
neutrals. The electron parallel momentum Eq. (6) now reads 



The ion dynamics is the same as above, so we use Eq. (8). The dynamics of neutrals is 
completely described by 




(21) 




(22) 




(23) 




(24) 
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This is used in Eq. (23), with the momentum conservation condition that now reads 
mnnnVne = meTioi^en, yielding 

Vezi = , a = IH . (25) 

Here, e = meno/ (ninnn), and in view of Eq. (2) this is a small quantity for any plasma. 
For instance, for an electron-proton plasma in a hydrogen gas m„ = = rUp, it is of the 
order 10~^ or less. Equation (25) is used in the electron continuity Eq. (3) yielding 

rii uo^ + idhiD ^ v'l^ 



uq uj -\- iak^D ' 
In the case a;, cu* ^ ^iD, Eq. (26) can be written as 



D = -^. (26) 



no 



]_ _j_ en cii. • — ^ ^ 



u;2 CO kjD kjD 



We have 1 » 2eUen/k'^D ^ ^vlJit?^ the real part is therefore very close to unity, while 
the imaginary part is simply iiuj — uj^)/klD. Consequently, combining Eq. (27) with 
Eq. (8), the real part of the frequency appears described as earlier by Eq. (14), while the 
imaginary part becomes similar to Eq. (16), viz. 



2 i^Pk'^n'^ - 



(28) 



klDuj^,{l + klpj) + k',cr 

However, we stress again that Eq. (16) follows from a formally incorrect procedure. 

Compared to the previously discussed e — i collision case (20), it is seen that i) in 
both cases the threshold is caused by the ion sound response, however, ii) the instability 
threshold in Eq. (28) is shifted towards higher frequencies (because kyPs is usually less 
than unity). 

Although the two cases describe two physically different plasma environments, the 
explanation for case ii) should be as follows. The lower threshold frequency in the e — i 
case implies that electrons experience a larger amount of collisions with ions within a wave 
period, which is in fact necessary to compensate for the ion movement in the parallel 
direction (due to the parallel electric field). This is because moving ions (in the same 
direction as electrons) represent a less efficient barrier for electron parallel motion and, in 
order to have the necessary phase shift between the density and potential, the electrons 
should have more collisions for the instabihty to take place. On the other hand, in the e—n 
case, a higher frequency (equivalent to a smaller amount of collisions) for the instability 
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to develop is possible because neutrals are less movable in the parallel direction (they 
do not react to the parallel electric field), and therefore they represent a more effective 
barrier. 

VII. To conclude, the self-consistent inclusion of the momentum conservation in the ion 
and electron equations, which originates from the collisions between the two fiuids, yields 
a different instabihty threshold that, to the best of our knowledge, has not been discussed 

in the literature so far. The correct expressions (20) and (28) should be used for estimates 
of the growth rate and the instability threshold for the coUisional drift mode. 
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Figure caption: 

Fig. 1. The evolution of the ion velocity from Eqs. (21, 22) for an electron-ion plasma, 
normalized to i^ieVo/ {uei + I'ie), and for Veo = Vo, Vio = 0. 



12 



